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SPREAD  F  THEORIES  -  A  REVIEW 


INTRODUCTION 

Our  understanding  of  equatorial  spread  F  (ESF)  phenomena  has  increased 
dramatically  over  the  past  few  years;  notwi thstand ing  the  fact  that  it  was 
discovered  over  four  decades  ago  (BOOKER  and  WELLS,  1938)  as  diffuse  echoes 
on  ionograms.  Much  of  this  advancement  in  understanding  comes  from  a  com¬ 
bined  theoretical  and  experimental  approach  to  the  problem.  Advances  in 
radar  backscatter  measurements,  satellite  and  rocket  in  situ  measurements, 
and  theoretical  and  numerical  simulation  techniques  have  provided  a  clearer 
picture  of  the  fundamental  mechanisms  causing  ESF.  The  study  of  ESF  pheno¬ 
mena  has  been  both  complicated  and  enriched  by  the  fact  that  the  attendant 
magnetic  field  aligned  ionospheric  irregularities  span  some  5-6  orders  of 
magnitude  in  scale  sizes.  The  present  paper  will  deal  with  ESF  theory  and 
at  that  only  those  theories  using  plasma  mechanisms  as  a  basis.  Experimen¬ 
tal  reviews  of  the  subject  can  be  founji  in  other  papers  in  this  edition. 

-Before  we  proceed  let  us  examine  the  basic  equatorial  nighttime 
ionospheric  F  region  geometry ./.  'Figure  1  depicts  this  geometry  under  which 
ESF  occurs.  N(y)  represents  the  background  zero  order  electron  density  as  a 
function  of  altitude  (y).  Gravity,  g,  points  down,  the  ambient  geomagnetic 
field,  B,  is  horizontal  and  pointing  north  (z),  and  k  represents  a  horizon¬ 
tal  perturbation  vector  and  points  westward  (x).  The  maximum  in  the 
electron  density  defines  the  F  peak.  The  bottomside  of  the  profile  steepens 
at  night  due  to  chemical  recombination  effects  and  upward  ambient  L  X  B 
elec.trodynamic  forces.  The  E  region  has  been  reduced  by  chemical 
recombination.  i  le  might  expect  this  equatorial  ionospheric  geometry  to  be 
unstable  to  a  variety  of  plasma  instabilities  because  in  the  plasma  physics 

community  fig.  1  represents  the  classical  flute  mode  geometry. 
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Notwithstanding  the  simple  geometry  exhibited  in  fig.  1,  several  basic 
questions  have  driven  researchers  to  seek  solutions,  oftentimes  not  simple, 
to  the  mysteries  of  ESF  phenomena.  Some  of  these  basic  questions  are  as 
follows.  (1)  Why  are  plasma  density  fluctuations  observed  on  both  sides  of 
the  F  peak  when  linear  theory  of  the  representative  instabilities  predict 
only  bottomside  growth  (unstable)?  (2)  Why  does  the  onset  of  topside 
spread  F  lag  the  onset  of  bottomside  spread  F?  (3)  How  are  the  large  plasma 
depletions  (bubbles)  produced?  (4)  What  is  the  power  spectrum  for  the 
plasma  irregularities?  (5)  What  is  the  relation  between  the  observed  back- 
scatter  radar  short  wavelength  irregularities  and  the  longer  wavelength 
irregularities  observed  by  in  situ  measurements  and  scintillation  measure¬ 
ments  ? 

At  this  juncture,  a  general  brief  review  will  be  given  of  ESF  theories. 

This  will  be  done  from  an  historical  perspective,  but  first  we  will  deal 
with  linear  theories  and  then  nonlinear  theories.  The  idea  that  ESF  could 
be  initiated  on  the  bottomside  of  the  F  region  (see  fig.  1)  by  a  Rayleigh- 
Tavlor  instability  was  first  proposed  by  DUNGEY  (1956).  DAGG  (1957)  sug¬ 
gested  that  ESF  phenomena  was  caused  bv  irregularities  produced  first  in  the 
E  region  and  then  coupled,  by  the  high  conductivity  along  the  magnetic 
field,  up  to  the  F  region.  It  was  first  suggested  by  MARTYN  (1959)  that  ESF 
was  a  manifestation  of  the  E  X  B  gradient  drift  instability.  CALVERT  (1963) 
proposed  that  the  downward  motion  of  the  nighttime  neutral  atmosphere  was 
responsible  for  ESF.  However,  this  mechanism  was  essentially  equivalent  to 
the  E  X  B  gradient  drift  mechanism  because  of  the  role  played  by  the  rela¬ 
tive  drift  between  ions  and  neutrals  in  determining  the  growth  rate  of  the 
instability.  It  should  be  noted  that  all  of  the  preceding  linear 
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irregularities  on  the  bottomside  of  the  F  region.  Tlie.se  theories  were  dis¬ 
missed  as  inadequate  by  FARLEY  et  al .  (1970)  because  they  could  not  explain 
the  wide  variety  of  3m  coherent  radar  backscatter  measurements  made  at 
Jicamarca. 

In  an  effort  to  address  the  fact  that  irregularities  were  observed  on 
both  the  bottomside  and  topside  of  the  F  region,  the  collisional  Rayleigh- 
Tavlor  instability  with  field  line  averaging  was  proposed  (BALSLEY  et  al., 
1972;  HAERENDEL,  1974).  By  averaging  (integrating)  the  density  along  the 
magnetic  field,  the  total  electron  content  profile  becomes  steeper  on  the 
bottomside  of  the  F  region  and  its  peak  is  raised  in  altitude  with  respect 
to  the  local  electron  density  peak.  Consequently,  this  allows  the  linear 
mechanism  to  operate  to  slightly  higher  altitudes,  ~  100  km  greater  alti¬ 
tudes,  but  still  would  not  explain  the  occurrence  of  irregularities  above  th 
new  peak.  In  the  wavelength  regime  30m  -  100m,  HUDSON  and  KENNEL  (1975) 
pointed  out  the  importance  of  the  collisional  drift  mode  in  ESF.  This  mode 
had  the  advantage  that  it  could  occur  on  both  the  top  and  bottomside  of  the 
F  region,  but  still  would  not  explain  the  longer  wavelength  irregularities. 
Finite  larmor  radius  (FLR)  corrections  were  also  applied  to  the  collision¬ 
less  and  collisional  Ray leigh-Taylor  instability  in  their  work.  These  FLR 
effects  had  more  of  a  stabilizing  influence,  at  short  wavelengths,  on  the 
col  1 ision 1  ess  rather  than  the  collisional  regime.  CHIU  and  STRAUS  (1979) 
reexamined  the  linear  theory  of  the  Ravleigh-Taylor  instability  in  the 
nighttime  equatorial  ionosphere,  in  the  100  km  to  the  bottomside  F  region 
altitude  range.  They  included  (a)  the  eastward  drift  of  the  ionosphere 
caused  by  the  nighttime  polarization  electric,  field;  (b)  the  eastward  night¬ 
time  neutral  wind;  and  (c)  recombination  in  the  E  and  F  regions.  They 
found  that,  well  below  the  bottomside  F  region,  the  Ravleigh-Taylor  mode 
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could  be  unstable,  being  driven  by  an  eastward  neutral  wind  (with  the  plasma 
density  gradient  in  altitude)  rather  than  by  gravity.  This  investigation 
dealt  with  instability  only  below  the  F  peak. 

Because  of  the  inability  of  linear  theory  to  explain  most  ESF 
phenomena,  several  nonlinear  theories  have  been  set  forth  in  an  attempt  to 
explain  the  various  observations.  HAERENDEI,  (1972;  1974)  was  the  first  to 
suggest  that  the  range  of  wavelengths  (many  kilometers  down  to  centimeters) 
observed  during  ESF  was  due  to  a  hierarchy  of  (multi-step)  processes.  The 
picture  portrayed  was  as  follows:  (i)  the  collisional  Rayleigh-Taylor 
(R-T)  instability  with  horizontal  wavevectors  is  driven  by  gravity  and  the 
vertical  background,  zero  order  electron  density  gradient  scale  length  on  the 
bottomside  (see  fig.  1);  then  (ii)  the  E  X  _B  gradient  drift  instability  with 
vertical  wavevectors  arises  due  to  the  horizontal  density,  large  amplitude 
variations  (with  shorter  scale  lengths  than  the  original  zero  order  vertical 
density  gradient)  set  up  by  the  collisional  R-T  instability;  then  (iii)  the 
inertia  (collisionless)  dominated  R-T  instability  arises;  and  finally  (iv) 
kinetic  drift  waves  grow  off  these  irregularities  after  they  reach  large 
amplitude.  Because  each  succeeding  process  is  driven  by  smaller  and  smaller 
plasma  density  gradient  scale  lengths,  all  scale  size  irregularities  will 
appear  almost  simultaneously. 

In  the  same  spirit,  HUDSON,  et  al.  (1973)  suggested  that  the  very 
smallest  scale  ('!  10m)  irregularities,  i.e.,  those  observed  by  coherent  back- 
scatter  radar  techniques,  weie  due  to  a  two  step  process.  In  this  prescrip¬ 
tion  a  long  wavelength  collisional  drift  mode  grows  on  the  zero  order  verti¬ 
cal  density  gradient  to  finite  amplitudes,  generating  horizontal  density 
gradients.  These  horizontal  gradients  in  turn  generate  an  unstable  kinetic 
short  wavelength  drift  mode  at  wavelengths  less  than  the  ion  gvroradius. 
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CHATURVKDT  and  KAW  (1975a, b)  have  examined  the  formation  of  steady-state  ESF 
irregularities  in  the  collisionless  R-T  regime,  using  a  nonlinear  one- 
dimensional  set  of  equations  to  describe  the  instability.  CHATURVF.IJI  and  KAW 
(1976)  interpreted  the  k  2  (where  k  is  the  wavenumber  perpendicular  to  the 
geomagnetic  field  and  equal  to  2tt/A  and  A  is  the  wavelength)  measured  power 
spectrum  of  the  F.SF  plasma  density  irregularities  in  terms  of  a  two  step 
theory.  In  this  particular  theory  longer  wavelength  R-T  modes  couple  to 
kinetic  collisional  drift  waves  in  such  a  manner  that  the  mode  coupling 
results  in  the  observed  k~2  spectrum. 

SCANNAPIECO  and  OSSAKOW  (1976)  performed  the  first  numerical  simulation 
of  the  long  wavelength  collisional  R-T  instability  mechanism  in  the  nighttime 
equatorial  F  region  ionosphere.  Tire  simulation  was  two  dimensional  (x,y  see 
fig.  1)  in  the  plane  perpendicular  to  the  ambient  geomagnetic  field,  using 
only  gravity  and  the  vertical  plasma  density  gradient  (with  an  actual  iono¬ 
spheric  electron  density  profile)  as  drivers  of  the  instability.  The  simula¬ 
tion  results  showed  that  tire  collisional  R-T  instability  generated  irregu¬ 
larities  and  plasma  bubbles  (plasma  density  depletions)  on  the  bottomside  of 
the  nighttime  F  region  which  then  rose  beyond  the  F  peak  by  nonlinear 
polarization  induced  E  X  15  motion.  This  was  a  major  breakthrough  as  it 
represented  the  first  theoretical  result  to  explain  how  long  wavelength 
irregularities  could  appear  on  both  the  bottomside  and  topside  of  the 
equatorial  F  region.  The  bubble  phenomena  was  in  accord  with  many  observa¬ 
tions  (KELLEY  et  al.,  1976;  WOODMAN  and  LA  HOZ,  1976;  McCLURE  et  al . ,  1977). 
An  analytic  nonlinear  mode-mode  coupling  theory  for  the  coherent  development 
of  the  collisional  R-T  instability  was  performed  by  CHATURVEDI  and  OSSAKOW 
(1977).  In  this  theory  linearly  unstable  collisional  R-T  modes  in  the  x-y 
plane  (see  fig.  1)  saturate  nonlinearlv  hv  exciting  linearly  damped  vertical 
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modes,  i.e.,  these  vertical  modes  rise  nonl inear Iv  to  large  amplitudes.  In 
addition,  the  theory  suggested  that  the  vertical  modes  would  be  dominant  and 
result  in  a  k-2  power  spectrum.  HUDSON  (1978)  extended  the  results  of 
CHATL'RVF.Dl  and  OSSAKOW  (1977)  to  the  collisionless  R-T  regime  and  reached 
similar  conclusions. 

Given  the  observational  and  numerical  simulation  evidence  for  ESF 
hubbies,  some  theoretical  research  has  been  directed  towards  bubble  pheno¬ 
menology  in  the  ESF  ionosphere.  Analytical  models  for  the  rise  of  ESF  bub¬ 
bles  in  the  collisional  and  collisionless  R-T  regime,  in  analogy  with  fluid 
bubbles,  were  given  by  OTT  (1978).  The  result  of  his  study  predicts  the 
bubbles  to  be  cylindrical  (circular  cap  at  the  top)  in  the  two  dimensions 
perpendicular  to  B.  In  the  collision  dominated  R-T  regime,  for  a  100% 
depletion,  the  bubble  rise  velocity  V„  =  g/v.  ,  where  v.  is  the  ion-neutral 
collision  frequency  and  g  is  gravity.  In  the  collisionless  or  inertia  domi¬ 
nated  R-T  regime,  for  a  100%  depletion,  the  bubble  rise  velocity  turns  out 
v 

to  be  V  ~  (Rg)2/2,  where  R  is  the  radius  of  curvature  at  the  top  of  the 
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bubble  so  that  this  velocity  becomes  bubble  size  dependent.  OSSAKOW  and 

CHATURVEDI  (1978)  have  presented  two  dimensional  (in  plane  perpendicular  to 

B) ,  analytic  models  for  the  rise  of  ESF  bubbles  in  the  collisional  R-T 

regime,  within  the  framework  of  the  electrical  analogy  with  barium  clouds. 

Different  bubble  geometries  were  investigated  with  the  results  of  all  studios 

expressible  as  V  =  (g/v.  )  f(6n/n  ),  where  f(6n/n  )  is  an  increasing  fune- 
Din  o  o 

tion  of  the  percentage  depletion  6n/n^  (with  the  case  of  a  100%  cylindrical 
depletion  being  g/v.  ).  SZUSZCZKWI CZ  (1  978)  has  discussed  the  chemistry  and 
transport  of  ESF  hubbies.  Investigation  of  vertical  induced  E  X  B  ESF  bubble 
motion,  incorporating  flux  tube  integrated  quantities  of  electron  content  and 
Pedersen  conductivity,  has  been  carried  out  by  ANDERSON  and  HAF.RF.NDEI,  (1979). 


The  si. lb  model  Lhat  thev  used  for  the  bubble  resulted  in  a  simple  expression 


for  tlu1  induced  electric  field  in  the  east-west  direction.  The  calculations 
showed  that  the  vertical  bubble  rise  velncitv  as  a  function  of  time  criti¬ 
cally  depended  on  the  background  electric  field.  Recent lv,  BURKE  (147b)  has 
investigated  the  physical  properties  of  ESE  bubbles  in  the  topside  1  region 
ionosphere  near  the  dawn  terminator.  The  studv  deals  with  the  effects  of 
the  K  region  after  sunrise  in  discharging  the  electric  fields  within  the  bub¬ 
bles.  The  results  show  that  the  conducting  K  region  can  effectively  halt  the 
upward  bubble  rise  velocity. 

COSTA  and  KELLEY  (L978a;b)  suggested  that  coherent  steepened  structures 
and  not  turbulence  would  give  a  k--  power  spectrum  for  ESE  irregularities. 

Furthermore,  these  sharp  plasma  density  gradients  could  cause  small  scale 
sizes  (~  20m)  by  collisionless  low  frequency  (<•<•  ion  gyrof  requency ,  ...) 
kinetic  drift  waves  through  a  two  step  process.  Tins  idea  was  set  forth  to 
explain  the  3m  Ticamarca  radar  observed  irregularities.  This  analysis  and 
subsequent  ones  (HiJRA  et  al . ,  1978;  HUBA  and  OSSAKOW,  1979a;b;  1980)  to 
explain  the  coherent  backseat  ter  radar  observed  irregularities  are  really 
linear  calculations  in  the  following  sense.  The  steepened  gradients  are 
formed  by  nonlinear  evolution  of  a  primary  instability,  e.g.,  R-T,  and  then 
one  performs  a  linear  theory,  for  a  secondary  instability  (e.g.,  kinetic 
drift  wave),  using  the  nonlinear  state  plasma  density  gradient  as  a  driver. 

Inclusion  of  el ec t ron- ion ,  electron-neutral,  and  ion-ion  collisions  by  HUBA 
and  OSSAKOW  (1979a)  into  the  low  frequency  drift  wave  theory  invoked  by  COSTA 
and  KEI.I.F.Y  (1978b)  showed  that  linear  theory  could  not  explain  the  3m 
irregularities.  Rather  one  could  get  down  to  only  about  10m  with  these  low 
frequence  drift  waves  (with  collisions). 

i : 
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HUBA  et  al.  (1978)  reported  on  very  small  scale  (lm  and  36  cm)  F.SF 


irregularities  which  were  experimentally  observed  by  coherent  radar  back- 
scatter  techniques  at  Kwajalein.  A  two  step  process,  utilizing  high 
frequency  O'-  ii^)  kinetic  drift  cyclotron  (DC)  or  lower  hybrid  drift  (LHD) 
instabilities  was  invoked  to  explain  these  irregularities  below  the  ion 
gvroradius.  This  linear  instability  analysis,  including  ion-ion  collisions, 
performed  on  a  nonlinear  state,  favored  the  formation  of  LHD  instabilities. 
The  results  provided  an  encouraging  explanation  of  the  very  short  wave¬ 
lengths.  More  detailed  analyses  including  various  collisional  effects  have 
been  performed  on  the  high  frequency  drift  waves  in  these  inhomogeneous 
plasmas.  HUBA  and  OSSAKOW  (1979b)  showed  how  cyclotron  resonances  could  be 
destroyed  in  weakly  collisional  inhomogeneous  plasmas.  Ion-ion  and  electron- 
electron  collisions  were  considered  and  the  results  showed  how  the  DC 
instability  went  into  the  LHD  instability.  SPERLING  and  GOLDMAN  (1980) 
looked  further  into  electron  collisions 3  effects  on  the  LHD  instability. 

They  originally  included  both  electron-ion  and  electron-neutral  collisions, 
but  in  their  final  analysis  only  showed  the  effects  of  electron-ion  colli¬ 
sions.  The  work  of  HUBA  and  OSSAKOW  (1980)  included  electron-neutral  as  well 
as  electron-ion  collisions  in  their  LHD  instability  explanation  of  11  cm 
irregularities  observed  at  high  altitudes  (~  500  km)  during  ESI'  conditions  at 
Kwajalein.  The  absence  of  these  irregularities  at  lower  altitudes  was  attri¬ 
buted  to  the  stabilizing  influence  of  electron-neutral  collisions.  All  of 
the  short  wavelength  (lm,  36  cm,  and  11  cm)  irregularities  observed  at 
Kwajalein  during  ESP  appear  to  be  explainable  (for  initiation  purposes)  by 
the  linear  theory  of  the  LHD  instability  in  the  condition  favoring  low  plasma 
densities  and  steep  plasma  density  gradients.  However,  the  region  encompas¬ 
sing  greater  than  one  meter  to  about  10 ra  seems  to  be  inaccessible  to  either 
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the  linear  theory  of  the  low  frequency  or  high  frequency  kinetic  instabili¬ 
ties  mentioned  above,  for  reasonable  plasma  density  gradients. 

KELLEY  and  OTT  (1978)  have  proposed  that  the  ESF  bubbles,  in  the 
eollisionless  R-T  regime,  generate  a  wake  with  vortices.  Two  dimensional 
fluid  turbulence  theory  was  then  applied  to  the  model.  This  resulted  in  the 
development  of  turbulence  at  shorter  and  longer  wavelengths  than  the  bubble 
size.  This  in  turn  led  to  a  prediction  of  k-1  for  the  power  spectrum  of  ESF 
irregularities  (which  does  not  appear  to  be  in  agreement  with  experimental 
observations  which  favors  k-2)  in  the  range  L^,-1  <  k  <  L^-',  where  is  the 
stirring  (bubble)  size  and  is  a  dissipation  length  cutoff. 

OSSAKOW  et  al .  (1979)  extended  the  earlier  numerical  simulation  work  of 
SCANNAPIECO  and  OSSAKOW  (1976)  on  nonlinear  ESF  to  study  the  dependence  on 
altitude  of  the  F  peak  and  bottomside  background  electron  density  gradient 
scale  length.  It  was  found  that  under  favorable  conditions,  i.e.,  high 
altitude  of  the  F  peak  and/or  steep  bottomside  background  electron  density 
gradients,  the  collisional  R-T  instability  caused  linear  growth  on  the  bot¬ 
tomside  of  the  F  region.  Subsequently,  this  caused  plasma  bubbles  (density 
depletions)  to  be  formed  on  the  bottomside  which  then  steepened  on  their  top 
and  rose  nonlinearly  above  the  F  peak  by  polarization  E  X  B  motion.  While 
rising  the  bubbles  elongated  in  the  altitude  direction.  This  sequence  of 
events  produced  irregularities  on  the  topside  of  the  F  region  where  linear 
theory  of  the  collisional  R-T  instability  would  predict  no  irregularities 
(stability).  High  altitude  of  the  F  peak,  small  bottomside  background 
electron  density  gradient  scale  lengths,  and  large  initial  bottomside  per¬ 
centage  depletions  yielded  large  vertical  bubble  rise  velocities.  Changing 
the  altitude  of  the  F  peak  (e.g.,  raising  it  80  km  above  the  case  of 
SCANNAPIECO  and  OSSAKOW  (1976))  exhibited  dramatic  effects  on  ESF  evolution. 
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The  higher  rise  velocities  due  to  the  elevated  F  peak  produced  results  in 
good  agreement  with  the  high  velocity  measurements  of  KELLEY  et  al .  (1976) 
and  MeCLURE  et  al.  (1977). 

Further  numerical  simulations  (7.ALESAK  and  OSSAKOW,  1980)  in  the 
collisional  R-T  regime  have  shown  that  horizontally  large  spatial  ESF  bubbles 
(~  50-100  km  wide),  with  almost  100%  depletions  (several  orders  of  magnitude 
change  in  plasma  density)  and  large  rise  velocities,  can  be  produced  with  an 
initially  long  wavelength  plasma  density  perturbation.  The  evolution  of  such 
bubbles  is  on  the  same  time  scale  as  the  small  scale  bubbles  studied  by 
OSSAKOW  et  al .  (1979),  but  the  plasma  comes  from  lower  altitudes  in  the 
equatorial  ionosphere.  This  last  fact  is  a  consequence  of  the  greater  pene¬ 
tration  of  the  polarization  fringe  field  into  the  lower  ionosphere  in  the 
large  versus  the  small  bubble  case  (see  next  theory  section).  KESK1NEN  et 
al.  (1980)  have  performed  nonlinear  numerical  simulations  of  the  intermediate 
wavelength  (100m  -  1  km)  collisional  R-T  instability  in  local  unstable 
regions  of  the  postsunset  bottomside  (~  300  km  altitude)  equatorial  F  region 
ionosphere.  It  is  found  that  the  linearly  unstable  modes  saturate  by  non¬ 
linear  excitation  of  linearly  damped  vertical  modes.  In  addition,  the  one 
dimensional  power  spectrum  for  these  irregularities  in  both  the  horizontal 
and  vertical  directions  is  «  k  n,  where  n  =  2-2.5.  These  results  support  the 
nonlinear  analytic  studies  of  CHATURVEDI  and  OSSAKOW  (1977).  Preliminary 
numerical  simulations  have  been  performed  (CHATURVEDI  et  al . ,  1978)  in  the 
collisional  R-T  regime  in  which  an  eastward  neutral  wind  (relative  to  the 
plasma  motion)  was  added.  These  results  indicated  that  the  wind  caused  a 
vertical  polarization  electric  field  which  in  turn  caused  the  ESF  bubbles  to 
move  westward  with  respect  to  the  plasma  in  agreement  with  the  observations 
of  MeCLURE  et  al.  (1977).  The  need  of  an  eastward  neutral  wind  to  produce 
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thi'  westward  drift  of  tin-  rising  bubbles  had  been  hypothesized  by  WOODMAN  and 


I -A  H07.  (197b),  OTT  (1978)  and  OSSAKOW  and  C.HATURVEDI  (1978). 

Much  work  and  significant  progress  in  the  theoretical  area  of  ESF 
phenomena  has  been  accomplished  in  the  past  few  vears.  In  spite  of  the  suc¬ 
cess  further  work  still  needs  to  be  done  (these  areas  are  summarized  in  the 
last  section  of  this  paper).  Theoretical  and  numerical  simulation  efforts 
are  continuing.  The  next  section  briefly  presents  some  of  the  representative 
theoretical  efforts  mentioned  in  the  preceding  paragraphs.  The  final  section 
presents  a  summary  concerning  F.SF  theory  and  remaining  problem  areas. 


THEORY 

In  this  section  we  present  some  representative  theoretical  and  numerical 
simulation  works  with  the  appropriate  references. 


Genera  1 


The  basic  plasma  fluid  equations  applicable  to  the  ESF  ionosphere  are  as 
foil ows : 
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12 


In  the  above  equations  the  subscript  u.  denotes  species  (e  is  electron,  i  is 
ion),  n  is  density,  V  is  velocity,  P  is  the  production,  v  is  the  chemical 
recombination  coefficient,  T  is  temperature,  V  is  the  gradient  operator,  B 
is  the  ambient  magnetic  field  (taken  to  be  uniform),  e  is  the  electronic 
charge,  c  is  the  speed  of  light,  m  is  mass,  £  is  gravity,  v.  is  the  ion- 
neutral  collision  frequency,  and  J  is  current.  In  addition,  we  have 
neglected  neutral  wind  effects.  Equation  (1)  is  the  continuity  equation, 

(2)  and  (3)  are  the  electron  and  ion  momentum  equations,  respectively,  (4) 
is  the  divergence  of  the  current  and  (5)  is  the  current  equation.  What  we 
have  in  mind  is  to  apply  the  set  of  equations  (1)  -  (5)  to  the  two  dimensions 
perpendicular  to  B  at  the  geomagnetic  equator,  making  various  approximations. 

If  inertial  terms  are  neglected  (i.e.,  the  left  hand  side  of  (3)),  we 
take  v  «  1  (valid  in  the  F  region),  and  we  make  the  quasi-neutral 

approximation,  ng  «  n.  &  n,  then  (1)  -  (3)  become  for  the  two  dimensions 
(x,y)  perpendicular  to  15 
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where  A.  =  eB/m.c,  B  =  Bz  and  we  have  set  T  =  T.  =  0  for  simplicity  (see 
ti—  e  i  1 

OSSAKOW  et  al . ,  1979).  We  have  neglected  the  £  x  z  motion  in  (7)  because 
compared  with  the  same  ion  motion  it  is  in  the  ratio  m^/mu  .  Making  the 
electrostatic  approximation  E  =  -  V<£  and  setting  <j>  =  $  +  <J> ,  where  t  lie  sub¬ 
script  zero  refers  to  equilibrium  or  zero  order  quantities  and  the  tilde 
refers  to  perturbed  quantities,  we  obtain 
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(10) 


In  arriving  at  (10)  we  have  put  (7)  and  (8)  into  (4)  and  (5)  and  obtained 

V 4>o  =  m.g/e.  Equation  (9)  is  the  electron  continuity  equation  in  the 

-  V ij>  x  B  drift  frame  (i.e.,  -  g  x  z/Q^  frame). 

To  obtain  a  linear  growth  rate  from  (9)  and  (10)  we  set  n  =  n  (v)  +  n 
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For  horizontal  waves,  i.e.,  k  =  k  this  reduces  to 
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Equations  (11)  and  (12)  represent  the  linear  growth  rate  in  the  collisional 
R-T  regime.  Linearizing  the  more  general  set  of  equations  (1)  -  (5),  with 
ion  inertial  terms,  would  have  resulted  in  the  growth  rate 

[-Vin  +  (Vin  +  48L)1^ 

Y  =  - __2 - Vr  03) 


which  reduces  to 
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Equation  (14a)  is  just  (12)  and  (14b)  represents  the  collisionless  (inertia 
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dominated)  R-T  regime. 


21)  Global  Numerical  Simulation  Results 


In  this;  section  we  will  outline  some  two  dimensional  (±B)  numerical 
simulation  results  of  the  collisional  R-T  instability  mechanism  (ZALESAK  and 
OSSAKOW,  1980).  The  simulation  plane  encompasses  an  altitude  extent  of  over 
200  km,  from  the  bottomside  to  topside  of  the  ESF  region.  Ionospheric 
profiles  of  n  (v),  v.  (v),  and  v  (y)  are  used  in  the  simulations  (see  figs. 

2  and  3).  A  simulation  mesh  of  140  (altitude,  y)  by  40  (east-west,  x)  was 
employed  with  Ay  =  2  km  and  Ax  =  5  km  (large  case,  L)  and  Ax  =  200m  (small 
case,  S) .  The  boundary  conditions  were  periodic  in  x  and  in  y  transmittive 
for  n  and  Neumann  for  the  induced  potential,  <J> .  Equations  (9)  and  (10)  were 
utilised,  with  n  obeying  dn  /dt  =  -  v  n  .  Two  different  tvpes  of  initial 
perturbations  and  different  wavelengths  corresponding  to  the  large  case  and 
small  case  were  used  (4  cases  in  all).  However,  here  we  will  only  present 
cases  IS  and  1L  which  are def ined  with  Ax  =  200m  and  Ax  =  5  km,  respectively, 
and  a  perturbation  of  the  form  (see  ZALESAK  and  OSSAKOW,  1980  for  more 
detail s) 
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In  contrast  to  this  more  localized  type  of  perturbation  which  lias  been  used 
before  (OSSAKOW  et  al . ,  1979)  the  pure  sinusoidal  one,  also  used  by  ZALESAK 
and  OSSAKOW  (1980),  will  not  be  shown  here. 
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RECOMBINATION 


e  ion-neutral  collision  frequency,  (solid),  and  recombination  coefficient  (rate),  vR  (dashed), 

function  of  altitude  used  in  the  numerical  simulations  (after  ZALESAK  and  OSSAKOW,  1980). 


Fig.  3  -  Sequence  of  four  plots  showing  isoelectron  density  contours  of 
calculation  IS  at  300,  700,  and  1000,  and  1200  sec.  Superimposed  on  each 
plot  is  a  dashed  line  showing  nQ  (y,t)  and  labeled  at  the  top.  Electron 
densities  are  given  in  cm  .  The  observer  is  looking  southward  so  that  B_ 
is  out  of  the  figure  (after  ZALESAK  and  0SSAKCW,  1980). 
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Figures  3  and  4  are  a  comparison  of  small  versus  large  horizontal 
extent.  Figure  2  shows  isodensity  contours  of  calculation  IS  at  times  300, 
700,  1000,  and  1200  seconds  after  initialization.  Figure  3  shows  the  same 
contours  at  the  same  times  but  for  calculation  If.  The  presence  of  much 
lower  density  plasma  fluid  in  the  bubble  in  calculation  1L  is  obvious.  There 
is  also  a  basic  difference  in  bubbie  morphology  at  late  times.  At  1200 
seconds,  IS  has  pinched  off  into  two  bubbles,  with  the  more  intense  one 
below  the  initial  central  bubble.  In  addition,  another  bubble  has  formed  in 
the  sides  of  the  computational  mesh.  These  structures  are  more  obvious  in 
the  plot  of  n(x,y)/nQ(y)  at  1200  seconds  for  IS  shown  in  fig.  5a.  The  maxi¬ 
mum  depletion  levels  are  70%  in  the  top  central  bubble,  97%  in  the  lower 
central  bubble,  and  95%  in  the  side  bubble.  In  the  lower  bubble  in  fig.  4a 
the  outermost  contour  has  n/no  =0.5  (50%  depletion)  and  the  innermost, 
fifth  contour  has  n/n^  =  0.03  (97%  depletion).  In  obtaining  percentage 
enhancements  and  depletions  one  subtracts  1.0  from  n/n^. 

Ttie  results  of  calculation  1L,  in  figs.  3  and  4,  show  a  single  plume  of 
depleted  ionization  at  1200  seconds,  with  no  secondary  central  bubble  or  side 
bubble.  In  fig.  5b  a  plot  of  n(x,y)/nQ(y)  for  case  1L  at  1200  seconds  is 
exhibited.  The  level  of  depletion  is  greater  than  99.9 %  for  the  entire 
10  km  by  70  km  oval  hole  located  inside  the  tenth  solid  contour  of  fig.  5b 
and  represents  at  least  a  three  ordersof  magnitude  decrease  (bite-out)  in 
plasma  density.  The  bubble  rise  velocity  for  each  case,  computed  from  the 
last  two  frames  of  figs.  2  and  3,  are  IS  =  210  m/s  and  1L  =  230  m/s.  These 
large  bite-outs  and  large  rise  velocities  are  in  agreement  with  the  observa¬ 
tions  of  McCLURF.  et  al  .  ,  (1  977). 

The  reason  that  the  bubble  is  more  depleted  in  the  11.  case  versus  the 
IS  case  (see  figs.  3,  4  and  5)  can  be  seen  with  the  aid  of  fig.  6,  which 
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Fig.  4 


Same  as  figure  3  but  for  calculation  lL  (after 
ZALFSAK  and  OSSAKOW,  1980) 
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value  of  n/nQ,  a  factor  of  0.5  times  the  previous  one.  The  first  enhancement  contour  is  for  n/n0 
=  2.0,  while  each  succeeding  contour  is  for  a  value  of  n/nQ,  a  factor  of  2.0  times  the  previous  one. 
The  superimposed  long  dashed  line  depicts  n  (y,t)  (after  ZALESAK  and  OSSAKOW,  1980). 
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zero  potential  is  not  plotted.  In  each  case  the  contours  are  chosen  to  space  the  twelve  contours 
(plus  the  zero  contour)  uniformly  from  maximum  to  minimum  potential.  The  large  dashed  line  depic 
n  (y)  at  700  s.  (after  ZALESAK  and  OSSAKOW,  1980). 


shows  tho  induced  electrostatic  potential,  !,  at  700  seconds.  Notice  that 
the  fringe  field  surrounding  the  bubble  in  the  II.  case  extends  to  lower 
altitudes  than  in  the  IS  case.  Since  the  contours  of  constant  potential  are 
in  fact  t low  streamlines,  more  plasma  from  lower  altitudes  is  being  drawn 
upward  in  the  11.  ease  than  in  the  IS  case.  The  effect  is  similar  to  the 
fringe  field  of  a  parallel  plate  capacitor.  In  that  case  the  distance  out¬ 
side  the  plates  that  the  electric  field  goes  is  proportional  to  the  plate 
separation.  Tn  the  case  depicted  in  fig.  6  the  role  of  plate  separation  is 
played  by  the  horizontal  size  of  the  bubble  (see  ZALESAK  and  0SSAK0W,  1980 
for  more  details). 

The  results  gleaned  from  these  numerical  simulations  of  large  horizontal 
wavelength  initial  perturbations  in  the  collisional  R-T  regime  are  as 
follows.  The  large  scale  length  initial  perturbations  evolve  non linearly 
into  large  horizontal  scale  length  ESF  bubbles.  These  bubbles  evolve 
approximately  on  the  same  time  scale  as  do  the  smaller  horizontal  scale 
length  counterparts.  The  plasma  comprising  these  large  bubbles  has  its 
origin  at  much  lower  altitudes  than  that  of  the  smaller  horizontal  scale 
lengths  bubbles  and  this  results  in  plasma  density  depletions  in  the  large 
hubbies  very  close  to  100%,  i.e.,  several  orders  of  magnitude  depletions, 
bocal  Numerical  Simulations 

Tn  this  section  we  will  show  an  example  of  some  computer  simulations  of 
intermediate  wavelength  (~  100m  -  1  km)  collisional  R-T  instability  evolution 
in  local  unstable  regions  of  the  postsunset  bottomside  (~  300  km  altitude) 
equatorial  F  region  Ionosphere  (KFSK1NEN  et  al.,  1980).  These  simulations 
are  called  local  compared  with  the  ones  presented  in  the  previous  section  in 
the  sense  that  a  chunk  of  the  F  region  is  looked  at  with  a  more  finely 


resolved  computational  mesh.  The  price  one  pays,  of  course,  is  that  one 


cannot  cover  a  great  altitude  or  east-west  extent  and  still  resolve  shorter 


wavelengths  (remember  the  F.SF  associated  irregularity  phenomena  cover  some 

5-6  orders  of  magnitude  in  scale  size).  In  these  particular  simulations  a 

64  by  64  computational  mesh  in  the  plane  perpendicular  (x,y)  to  Ji  is  utilized 

with  Ax  =  Ay  =  15m  (so  we  basically  have  a  1  km  by  1  km  chunk).  Simulations 

have  been  run  for  ambient  bottomside  electron  density  gradient  scale  lengths 

I.  =  5,  10,  and  15  km  (where  we  have  taken  n  (y)  =  N  (1  +  (v/L)).  Bottomside 

ionospheric  parameters  corresponding  to  an  altitude  of  300  km  has  been  used. 

With  the  above  ideas  in  mind,  equations  (9)  and  (10)  have  been  used, 

except  that  a  term  v^n^  has  been  added  to  the  right  side  of  (9).  If  one 

sets  n  =  n  in  (9)  this  shows  that  dn  /3t  =  0,  i.e.,  n  is  taken  not  to 
o  o  o 

change.  In  addition,  the  simulations  were  initialized  with  a  two  dimensional 
perturbation  of  the  following  form 

n  (x. > >  ’  ? - 5}-  =  A  sink  x  cos  k  y  +  A,  nsin  2  k  y  (16) 

nQ  1,1  x  yJ  2,0  yJ 

which  is  the  form  suggested  in  the  analytic  nonlinear  mode  coupling  work 

of  CHATURVEDI  and  OSSAKOW  (1977).  In  addition,  the  following  values  were 

taken:  k  =  k  2ir/960m,  A,  ,  =  10-4,  and  A„  .  =  2  x  10-e. 
v  x  l ,  1  2,0 

The  nonlinear  aspects  of  these  simulations  for  the  L  =  10  km  case  are 
shown  in  figs.  7-9.  Figure  7a  shows  an  isodensity  contour  plot  of  n/n^  at 
t  =  0.  The  contours  describe  a  sequence  of  enhancements  and  depletions 
arranged  in  a  checkerboard  fashion  in  line  with  the  initial  perturbation 
given  by  (16).  Figures  7b  and  7c  show  the  time  evolution  of  the  initial 
perturbation  at  t  =  3000  and  6000  sec. ,  respectively .  The  density  fluctuation 
contours  at  t  =  10,000  sec.  are  depicted  in  fig.  7d  and  here  some  vertical 
elongation  and  steepening  can  be  observed  along  with  smell  scale  structures. 
In  figs.  8a  and  8b  representative  one  dimensional  horizontal  P(k^)  and 


23 


r 


(a) 


(b) 


•J 

k  / 


I, 


i.. i.  .  ‘.j..  .. .. ,  ii.«. .  V.41:  .  ..i  .  .1. 


(c)  (d) 

Fig.  7  -  Isodensity  contours  of  n/n0  for  L  =  10  km  at  (a)  t  =  0,  (b)  t  =  3000  s., 
(c)  6000  s.,  and  (d)  t  =  10,000  s.  Solid  contours  denote  n/nD  >0;  dashed  con¬ 
tours  denote  n/nD  <0.  All  contours  are  evenly  spaced  with  the  y  axis  vertical 
and  the  x  axis  horizontal.  The  tick  marks  represent  grid  point  location  (after 
KESKINEN  et  al.,  1980). 
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Fig.  8  -  One  dimensional  (a)  horizontal  PO^)  and  (b)  vertical  P(ky)  power 
spectra  versus  kx  and  ky,  respectively,  for  1  *  10  km  at  t  =  10,000  s.  The 
solid  curves  are  the  least  squares  fit  to  the  results  from  the  numerical  sim 
ulations  (dots),  and  k_  =  2*/ 960m  is  the  fundamental  wave  number. 


vertical  P(k  )  power  spectra,  in  the  nonlinear  late  time  regime 

(t  =  10,000  sec.),  have  been  plotted.  In  'his  saturated  regime  we  see  that 

P(k  )  “■  k  and  I’(k  )  ^  k  for  A  ~  70  -  960m.  These  spectral 

x  x  v  y 

indices  are  in  agreement  with  experimental  observations  (DYSON  et  al . ,  1974; 

COSTA  and  KELLEY,  1978a),  global  numerical  simulations  (SCANNAP1EC0  and 

OSSAKOW,  1976)  and  nonlinear  analytic  theory  (CHATURVED1  and  OSSAKOW,  1977). 

Figure  9  depicts  the  evolution  of  the  A  and  A  modes.  Initially,  the 

1,12,0 

A,,  q  vertical  mode  is  damped  as  predicted  by  linear  theory.  As  the  linearly 
unstable  A  mode  increases  it  nonlinearly  triggers  the  growth  of  A7  .  At 

I  ,  i  /  ,  u 

late  times,  these  two  modes  arrange  themselves  in  accord  with  the  nonlinear 

theoretical  predictions  of  CHATURVEDI  and  OSSAKOW  (1977).  From  fig.  9,  we 

note  that  in  the  quasi-final  state  A  „  <  A,  , ;  wnereas,  for  the  L  =  5  km 

2,01,1 

case  depicted  in  KESKINF.N  et  al .  (1980)  A  >  A  (the  corresponding 

2,0  1,1 

figures  in  KESKINEN  et  al . ,  1980,  to  those  shown  here  as  figs.  8  and  9,  are 
for  L  =  15  and  L  =  5  km,  respectively  and  so  are  not  the  L  =  10  km  case). 

CHATURVEDI  and  OSSAKOW  (1977)  showed  that  for  a  two  dimensional 
perturbation  of  the  form  (16),  the  dominant  nonlinearity  in  the  continuity 
equation  is  the  XgxB  *  term’  Taking  n/n^  =  Aj  ^  sin(kxx-uit)cos  k  y  + 

A„  n  sin  2k  y  gives  the  equation  for  A  and  A 

■*cjU  V  *  ,  1  /-,U 


dA 


1,1 


at 


Yi,i  Ai,i 


2  a  A  A 

1,1  2,0 


(17) 


dA 

,0 


at 


(18) 


where 


Tl 


1 


2v .  L 


i  n 


n 

o 


an 

_ o 

ay 


(19) 
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(20) 


k2  =  k2  +  k2 
x  y 


(21) 


(22) 


(23) 


It  is  (22)  and  (23)  which  are  given  as  Q  and  in  fig.  9  (with  Vin  - 

0.6  s-1 ,  v  =  3  x  10~4  s_1) . 

K 

The  results  gleaned  from  these  intermediate  wavelength  collisional  R-T 
instability  numerical  simulations  in  the  nighttime  equatorial  F  region  are: 
(i)  linearly  unstable  modes  saturate  by  nonlinear  excitation  of  linearly 
damped  vertical  modes;  (ii)  a  vertical  elongation  is  observed  along  with 
small  scale  structure;  and  (iii)  the  one  dimensional  vertical  and  horizontal 
power  spectra  for  the  plasma  density  fluctuations  a  k  n ,  where  n  =  2-2.5. 

FSF  Simulation  Results  in  General 

We  can  summarize  the  nonlinear  numerical  simulation  results,  to  date, 
as  follows.  The  collisional  R-T  instability  can  cause  linear  growth  on  the 
bottomside  ot  the  ESF  region.  Plasma  depletions  (bubbles)  form  on  the  bot- 
tomside,  steepen  on  their  top  and  nonlinearlv  rise  to  the  topside  above  the 
F  peak  by  polarization  E  X  B  motion.  High  altitude  of  the  F  peak,  small 
bottomside  background  electron  density  gradient  scale  lengths,  and  large  per 
rentage  depletions  yield  large  vertical  bubble  rise  velocities,  with  the 
first  two  conditions  favoring  collisional  R-T  growth.  Horizontally  large 
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spatial  bubbles  (~  50  -  100  km  widths)  can  bo  produced  by  long  wavelength 
perturbations.  Preliminary  results  (CHATURVEDI  et  al.,  1978)  show  that 
westward  drifts  of  the  bubbles,  with  respect  to  the  main  plasma,  can  be  pro¬ 
duced  by  an  eastward  neutral  wind.  Results  of  intermediate  wavelength 
(100m  -  1  km)  simulations,  of  the  collisional  R-T  instability  in  local 
regions,  show  saturation  by  nonlinear  generation  of  linearly  damped  vertical 
modes  and  a  one  dimensional  power  spectrum  “  k  n,  where  n  =  2-2.5. 

Very  Small  Scale  10m)  Irregularities 

We  will  not  dwell  much  in  this  area  because  another  paper  in  this, 
edition  can  be  found  on  the  subject.  Some  basic  ideas  will  be  presented. 

The  work  in  this  section  is  being  applied  to  account  for  the  coherent  back- 
scatter  radar  measurements  made  of  small  scale  irregularities  at  various 
equatorial  locations  (3m  at  Jicamarca,  Peru;  lm,  36  cm,  and  11  cm  at 
Kwajalein,  Marshall  Islands).  The  calculations  are  multilinear,  using 
kinetic  theory  and  result  in  various  plasma  kinetic  drift  wave  modes.  Kinetic 
theory  is  used  because  of  the  small  sizes  involved  (remembering  that  the  ion 
gyroradius  is  ~  6m  for  the  ESF  ionosphere).  It  is  multilinear  because  it 
depends  on  a  two-step  process  whereby  linear  theory  is  performed  on  the  non¬ 
linear  s-r.-ite.  By  using  linear  theory  we  are  attempting  to  account  for  the 
radar  observed  irregularities  by  the  simplest  means,  i.e.,  linear  excitation. 
This  does  not  preclude  their  being  excited  by  truly  nonlinear  means.  How¬ 
ever,  if  they  can  be  excited  in  a  linear  fashion  this  will  give  us  a  warm 
feeling.  The  driving  density  gradients,  in  these  calculations,  are  thought 
to  arise  through  a  primary  instability  (e.g.,  long  wavelength  collisional 
R-T  instability),  driven  by  the  zero  order  background  ionospheric  equatorial 
F  region  vertical  electron  density  gradient,  achieving  a  large  amplitude 
state.  The  zero  order  ionospheric  electron  density  gradient  is  of  larger 
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scale  size  than  the  primary  instability  electron  density  gradient  scale 
size,  which  would  be  of  the  order  of  the  instability  wavelength.  Because  the 
calculations  are  kinetic  they  employ  particle  distribution  functions,  which 
would  be  described  by  a  Vlasov  equation  in  the  col  1  is i on  1  ess  limit  (COSTA  and 
KELLEY,  1978b)  and  contain  an  additional  collision  operator,  i.e.,  a  Fokker- 
Planck  type  equation,  in  the  collisional  regime  (see  HUBA  and  OSSAKOW, 
1979a;b;  1980;  SPERLINC  and  GOLDMAN,  1980). 

We  will  present  some  general  concepts  regarding  kinetic  drift  waves 
which  will  be  useful  for  the  calculations  which  have  been  done.  The  basic 
geometry  is  such  that 

B  =  B  z  (24) 


no  =  no(x)  (25) 

(similar  to  the  zero  order  background  equatorial  ionospheric  geometry)  where 
here  nQ(x)  arises  due  to  the  primary  instability.  The  orthogonal  coordinate 
system  is  completed  by  the  ion  and  electron  diamagnetic  velocities  being 
along  the  y  axis. 


V,  =  (V,.  -  V,  )  y 
^-d  di  de  J 

V , .  =  (V.2/2fi.)  (Sin  n  /Sx) 
d  l  it  o 


(26) 

(27) 


V ,  =  -  (V  2/ 2C5  )  (Sin  n  /Sx)  (28) 

de  e  e  o 

i„ 

where  V.  =  (2T.  /m.  )2  and  the  larmor  radius  is  defined  by  r  =  V /Q  . 

i  ,e  i  ,e  i  ,e  a  a  a 

The  linear  analysis  is  then  performed  with  this  in  mind. 

First  the  low  frequency,  uu  «  11 . ,  collisionless  drift  wave  calculations 
of  COSTA  and  KELLEY  (1978a; b)  for  k  r^  ~  1  were  performed  to  explain  the  3m 
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ESF  irregularities  observed  by  the  Jicamarca  radar  (50  MHz)  in  the  coherent 
backseat  ter  mode.  HUBA  and  OSSAKOW  (1979a)  included  electron-ion,  electron- 
neutral,  and  ion-ion  collisions  in  the  low  frequency  drift  waves  and  showed 
that  linear  theory  could  not  explain  the  3m  regime.  The  theory  could  only 
show  instability  down  to  ~  10m  for  reasonable  ESF  plasma  densities  and 
density  gradient  scale  lengths. 

In  a  series  of  papers,  the  linear  theory  of  high  frequency,  ul>  Q ^ , 
drif  l  waves  (HUBA  et  al . ,  1978;  HUBA  and  OSSAKOW,  1979b;  1980;  SPERLING  and 
GOLDMAN,  1980)  for  kr..  ^  1  and  kr^  ~  1  have  been  investigated  including,  at 
various  stages,  ion-ion,  electron-electron,  electron-ion,  and  electron- 
neutral  collisional  effects.  These  studies  have  been  undertaken  to  account 
for  the  lm,  36  cm,  and  11  cm  ESF  irregularities  observed  at  Kwajalein  with 
the  ALTAIR  (155  MHz  and  415  MHz)  and  TRADEX  (1320  MHz)  radars  operating  in 
the  coherent  backscatter  modes.  These  wavelengths  appear  to  be  explainable 
by  the  linear  theory  of  the  high  frequency  lower  hybrid  drift  instability 
favoring  low  densities  and  steep  plasma  density  gradients.  Consequently, 
within  the  context  of  linear  theory  while  these  high  frequency  instabilities 
can  go  up  to  ~  lm  and  the  low  frequency  instabilities  can  go  down  to  ~  10m, 
there  is  a  gray  area  between  l-10m  that  is  inaccessible  by  linear  theory. 

CONCLUSIONS  AND  SUMMARY 

The  basic  physical  picture  of  ESF  that  emerges  is  as  follows.  After 
sunset  the  E  region  begins  to  recombine  and  there  is  effectively  no  E  region 
conductivity  to  short  out  any  irregularities  in  the  F  region.  Due  to  recom¬ 
bination  and  possible  electrodynamic  effects  (upward  motion  of  the  F  region 
plasma  due  to  E  X  B)  the  bottomside  F  region  background  electron  density 
gradient  begins  to  steepen.  The  steepening  due  to  plasma  E  X  B  uplift 
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caused  by  an  eastward  electric  field  is  much  like  barium  cloud  steepening  in 
equatorial  geometry  (note  a  downward  motion  of  the  neutral  atmosphere  will 
have  the  same  effect).  When  the  altitude  of  the  F  region  is  high  enough  and/ 
or  the  bottomside  background  electron  density  gradient  steep  enough  to  over¬ 
come  recombination  effects,  plasma  density  fluctuations  will  begin  to  grow  on 
the  bottomside  via  the  collisional  Rayleigh-Taylor  instability  mechanism  (or 
possibly  initiated  by  the  E  X  B  gradient  drift  instability;  in  any  case  a 
fluid  type  gradient  instability).  These  irregularities  will  in  turn  form 
plasma  density  depletions  (bubbles)  on  the  bottomside  which  will  then  non- 
linearly  rise  by  polarization  E  X  13  motion  through  the  F  peak  and  cause 
topside  spread  F.  (The  higher  the  F  peak  and/or  the  sharper  the  bottomside 
gradient,  and  the  higher  the  %  depletion  the  faster  the  rise  rate  of  the  bub¬ 
ble  and  the  faster  the  evolution  of  ESF.)  These  steepening  bubbles  (on  their 
topside)  can  then  bifurcate  and  form  smaller  scale  structure  by  a  cascade  or 
two-step  mechanism.  An  eastward  neutral  wind  can  cause  the  bubbles  to  move 
westward  with  respect  to  the  bulk  plasma  motion.  The  long  wavelength  irregu¬ 
larities  have  a  power  law  PSD  ^  k~2i? .  It  is  presently  believed  that  the 
small  scale  (^  3m)  radar  backscatter  observed  irregularities  (kinetic  tvpe 
instabilities)  are  created  by  the  steep  plasma  density  gradients  created  by 
the  primary  long  wavelength  fluid  type  instabilities  (collisional  Rayleigh- 
Taylor,  etc.).  The  turn  on  and  turn  off  of  these  very  small  scale  irregu¬ 
larities  and  their  connection  to  the  longer  wavelength  irregularities  is  not 
well  understood  at  present  and  requires  further  investigation. 

Although  much  work  has  been  done  in  recent  years  on  ESF  and  a  coherent 
picture  is  emerging,  there  are  still  remaining  problems  worthy  of  investiga¬ 
tion.  Some  of  these  problems  are  as  follows.  (1)  The  relation  of  the  very 
small  scale  radar  observed  irregularities  to  the  longer  wavelength 
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irregularities  is  not  well  understood.  (a)  Quantitatively,  what  are  the 


turn  on  and  turn  off  criteria  for  these  very  small  scale  lengths?  (b)  A  more 
complete  study  of  the  collisional  effects  on  plasma  kinetic  drift  wave  insta¬ 
bilities  is  required.  (c)  A  determination  of  the  nonlinear  saturation  of  the 
small  scale  irregularities,  including  saturated  amplitudes  and  power  spectra, 
is  needed.  (2)  The  effects  of  ion  inertia,  including  the  investigation  of 
topside  structure,  plasma  bubble  evolution  into  collisionless  altitudes 
(v^  «  4g/L) ,  nonlinear  saturation,  and  bubble  pinch  off  or  wedge  formation 

has  to  be  studied.  (3)  We  need  to  do  an  investigation  of  the  effects  of 
neutral  winds,  electric  fields,  and  changing  initial  conditions  on  our  pre¬ 
vious  plasma  fluid  long  wavelength  studies  (particularly  the  numerical 
simulations).  (4)  An  analytic  description  of  many  ESF  bubbles  is  lacking. 
Also,  how  do  they  decay  and  what  is  the  role  of  diffusion?  (5)  There  is  a 

need  to  study  three  dimensional  effects,  including  the  effects  of  kji  (wave- 

>1 

number  parallel  to  B ) .  (6)  The  effects  of  other  regions  of  the  ionosphere, 

e.g.,  the  E  region,  needs  to  be  investigated.  These  problem  areas  should  keep 
us  busy  for  the  foreseeable  future. 
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